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Abstract 

One of the leading issues in high-Tc copper oxide superconductors is the origin of the pseudogap phase in the 
underdoped regime of their phase diagram. Using polarized neutron diffraction, a novel magnetic order has been 
identified as an hidden order parameter of the pseudogap as the transition temperature corresponds to what is 
expected for the pseudogap. The observed magnetic order preserves translational symmetry as predicted for orbital 
moments in the circulating current theory. Being now reported in three different cuprates families, it appears as a 
universal phenomenon whatever the crystal structure (with single Cu02 layer or bilayer per unit cell). To date, it 
is the first direct evidence of an hidden order parameter characterizing the pseudogap phase of high-Tc cuprates. 
To cite this article: P. Bourges and Y. Sidis, C. R. Physique xx (2010). 

Resume 

Ordre magnetique de Petat pseudogap des oxydes de cuivre supraconducteurs a haute temperature 
critique Un des problemes majeurs des oxydes de cuivre supraconducteurs a haute temperature critique est 
I'origine de la phase pseudogap dans I'etat sous-dope. En utilisant des mesures de diffraction de neutrons polarises, 
nous avons identifie un ordre magnetique cache qui apparait a la temperature attendue de cet etat de pseudogap. 
Get ordre magnetique respecte la symetrie de translation du reseau comme cela a ete predit dans la theorie des 
boucles de courants circulants. Nous avons generalise cette decouverte dans trois families distinctes de cuprates. 
Cette mesure est la premiere preuve experimentale directe d'un parametre d'ordre universel de I'etat de pseudogap 
des oxydes de cuivre supraconducteurs. 
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1. Introduction 

The origin of high-Tc superconductivity in copper oxide materials is still hotly debated more than twenty years 
after its discovery. On the one hand, conventional phonon-mediated superconductivity has been advocated to 
explain anomalies in electronic spectroscopies although the electron-phonon coupling seems insufficient to explain 
the high value of the critical temperature. On the other hand, the antiferromagnetic (AF) spin fluctuations 
observed in these strongly electronic correlated materials could also lead to an unconventional superconducting 
(SC) pairing mechanism, where the conventional electron-phonon coupling would be replaced by a spin-fermion 
coupling. For instance, unconventional AF excitations, the so-called resonance modes, have been reported in the 
SC state of most of cuprates families [1]. Assuming a rather large spin-fermion coupling (~1 eV), these collective 
AF excitations may account for several anomalies in charge excitation spectra and for the angular momentum 
dependence of the SC gap (see for instance Ref. [2]). However, the spin fluctuation mediated pairing scenario has 
to face a serious problem: high-Tc superconductivity survives even when the AF spin fluctuation spectral weight 
becomes strongly reduced. If certain aspects of the physics of cuprates can be understood using either phonons 
or AF spin fluctuations, the mystery of high-Tc superconductivity in cuprates remains unsolved. 

All the high-Tc SC cuprates share a common crystallographic structure: They are layered materials characterized 
by the stacking of Cu02 planes. These planes are paved with squared Cu02 plaquettes. The charge density in 
Cu02 planes can be tuned using either electron or hole doping. All the hole doped high-Tc SC cuprates exhibit 
the same remarkable phase diagram (Fig. 1). These compounds are AF Mott insulators at zero doping. The AF 
state is quickly destroyed once a small amount of doped holes is introduced in the Cu02 planes. Increasing the 
hole doping, the system becomes metallic and superconducting below the critical temperature Tc. At optimal 
doping (p~0.16 holes /Cu), Tc reaches its maximum value. Two distinct regimes develop on both sides of the 
optimal doping. In the overdoped regime when increasing the hole doping, the electronic properties in the normal 
state can be described using a moderately correlated Fermi liquid picture. At variance, in the underdoped regime 
when reducing the hole doping, the materials behave as a strongly correlated metals and the standard Fermi 
liquid picture fails to account for their unconventional electronic properties. In particular, the underdoped regime 
is dominated by a phase with highly unusual physical properties where magnetic, transport and thermodynamic 
measurements point towards a diminution of the electronic density of states below a temperature T* [3,4,5] 
although the system remains metallic. As the density of states remains non-zero at the lowest temperature, this 
phase has been named as a pseudogap phase since its first evidence[6]. However, below T*, a q-dependent gap opens 
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Figure 1. Generic phase diagram of high temperature superconducting copper oxides as a function of hole doping (p). In the Mott 
insulation state, the antiferromagnetic phase develops below the Neel temperature (Tn). Upon doping, the system becomes metallic 
and superconducting (SC) below the superconducting critical temperature (Tc). In the normal state, the lost of spectral weight on 
certain portions of the Fermi surface highlights the opening of the mysterious pseudogap phase below the temperature T*. 
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in the single-particle excitation spectrum measured in angle-resolved photoemission (ARPES) experiments with a 
maximum at the wave vector M=(7r,0) and symmetry related points of the Brillouin zone. Other spectroscopies [4] 
spanning from scanning tunneling microscopy (STM), electronic Raman spectroscopy and optical conductivity 
reported anomalies at similar energy. The understanding of the microscopic origin of the pseudogap phase and of 
its interplay with unconventional d-wave superconductivity has become one of the major challenges to overcome 
in order to crack the mystery of superconducting cuprates. 

Two major classes of theoretical models attempt to explain the pseudogap state. It has been at first proposed 
that the pseudogap phenomenon can be described as a precursor of the superconducting d-wave gap, but with no 
phase coherence, which would occur only below Tc [7]. This preformed pairs scenario has largely been favoured 
over the years in relation to the doping dependence of the SC gap which does not follow Tc (as it is expected 
in the conventional BCS theory of superconductivity) but rather T*. This has been discussed in conjunction of 
many experimental techniques such as ARPES, STM, electronic Raman scattering experiments. However, due to 
lacking experimental evidence of such preformed pairs [8], this approach is questioned. 

Many other theories attribute the pseudogap origin to the proximity of a competing state, but there is a 
wide disagreement about the nature of this state. In certain theoretical models, the pseudogap state is a long 
range ordered phase, with a well defined order parameter and a related broken symmetry [9,10,11,12]. In this 
scenario, the ordering temperature T* and the order parameter should vanish at a quantum critical point (QCP) 
located behind the SC dome, close to optimal doping (see dashed line in Fig. 1). A large number of experimental 
properties ranging between transport, thermodynamic and magnetic [5] point towards the existence of a QCP 
for a doping level of p=0.19. The order parameter may involve charge and spin density waves or charge currents 
fiowing around or inside Cu02 plaquettes [9,10,11,12]. Interestingly, while the low-doping phase may compete with 
superconductivity, the fluctuations associated with the broken symmetry are expected to be rather strong around 
the QCP and could play the role of pairing glue leading to high temperature superconductivity. This appealing 
scenario has nevertheless to face two major experimental facts: there is no clear jump of the specific heat at T*, 
as expected for many phase transitions (but not all) and there is no indication that the translation symmetry of 
the lattice (TSL) is broken in the pseudogap state. 

Alternatively, other theoretical approaches have been developed still on the basis of a competing electronic 
instability resulting from strong electronic interaction, but in these approaches the pseudogap state is no longer an 
ordered phase [13,14,15,16,17]. It is a disordered phase dominated by the fiuctuations associated with a competing 
electronic instability and T* is only a crossover of dynamical properties. Nevertheless, a true ordered state could 
be stabilized at low temperature at the expense of superconductivity when applying an external perturbation 
(external magnetic field, disorder, structural distortion, etc.). Such an approach is well illustrated by the stripes 
model where charges self-organize to form lines, separated AF domains in anti-phase [13]. When ID charge stripes 
appear, but still fluctuate, the C4 rotational symmetry is flrst broken. When fluctuating stripes are pinned down 
on the lattice or defects, they become static and the resulting charge and spin order ultimately breaks the TSL. In 
La2-xSrxCu04 (LSCO), incommensurate spin fluctuations are usually interpreted in terms of fluctuating stripes, 
that can be pinned down at low temperature by disorder or stabilized under an external magnetic field [18]. A 
static stripe phase has been further reported in La2-xBaxCu04 or in (La, Nd)2-xSrxCu04 around the critical 
hole doping p=l/8. In addition, in strongly underdoped YBa2Cu3 06.45 (p=0.085) [19] as well as for similar 
dopings [20], anisotropic incommensurate spin excitations are observed at low temperature, suggesting that the 
C4 rotational invariance of the system is spontaneously broken below ^ 150 K. These fluctuations further freeze at 
very low temperature, yielding a glassy short range spin density wave (SDW) state. This SDW state is enhanced by 
applying a weak external magnetic field [21], but still remains at short range. Meanwhile, small Fermi pockets have 
been detected in quantum oscillation measurements carried out at high magnetic field [22,23,24] in YBa2Cu306+x 
system (only for x>0.5, i.e, p>0.09). The most simple explanation for these Fermi pockets is a Fermi surface 
reconstruction[23,25], triggered by a spin and/or charge density wave order breaking the TSL. The reconstruction 
of the Fermi surface inferred from quantum oscillation measurements could be linked to the SDW ordering reported 
in YBCO with and without magnetic field [19,21,20], but for a lower hole doping (x<0.5, i.e, p<0.09). However, 
one should be rather careful since the correlation length associated with the glassy SDW state are too short to 
account for the observed quantum oscillation and no quantum oscillation have been reported so far at the small 
hole doping at which the SDW is observed. 

In following sections, we will present compelling evidence of an ordered magnetic phase in the pseudogap state. 
We will show that this order phase displays strong similarities with the circulating current (CC) phase proposed 
by CM. Varma [10,11]. The CC-phase is associated with a Q=0 electronic instability and therefore preserves the 
TSL. Furthermore, its generalized Ising (Ashkin- Teller) phase transition does not produce any strong jump in the 
specific heat for a large range of parameters. While a 3D long range Q=0 magnetic order is clearly visible in the 



3 



pseudogap state of system such YBa2Cu306+x (YBCO) or HgBa2Cu04+(5 (Hgl201), we will show that this order 
is frustrated in LSCO and is likely to be in competition with a stripe instability. 



2. Evidencing circulating currents 

2.1. Loop current order 

Beyond usual charge or spin instabilities, more exotic electronic phases, spanning various CC states [9,10,11,12], 
have been proposed to account for a hidden order parameter associated with the pseudogap. One may have a single 
charge current per Cu02 plaquette staggered in the neighbouring cell, referred to as a D-density wave (DDW)[12]. 
The DDW state implies a doubling of the unit cell. The interest for such a state has been recently revived [26,27] 
owing to the observation of small Fermi pockets in quantum oscillation measurements at high field. In the DDW 
state, charge loops give rise to orbital-like magnetic moment perpendicular to Cu02 planes. The AF arrangement 
of such orbital-like moments should be directly probed by neutron scattering diffraction. However, to date, the 
neutron detection a DDW state is rather scarce [28]. 

Other CC states have been predicted [9,10,11] from the 3-band Hubbard model involving both copper d-orbitals 
and in-plane oxygen p-orbitals. As a result of the Cu-0 repulsion term, a loop-current state is stabilized yielding 
patterns with two (CC — Bn phase), or four (CC — Bi phase) current loops per unit cell. The phase CC — Bn 
corresponding to two opposite current loops is depicted in Fig. 2. A, it belongs to the Eu irreductible representation 
of the D4h point group [30,31]. A current flows from the Cu atom through the nearest oxygen atoms, then back 
to Cu. By principle, these phases break the time-reversal symmetry (TRS). The broken TRS has been indeed 
observed by ARPES in Bi2Sr2CaCu2 08+5 (Bi2212) through a spontaneous dichroism below T* [29]. As sketched 
in Fig. 2.B, each closed current loop produces an orbital magnetic moment which should be measurable by neutron 
diffraction. A sizeable orbital magnetic moment of 0.1 jiB is typically expected [9], pointing perpendicularly to 
the Cu02 planes. In contrast to the DDW state, these phases preserve the TSL, as the same pattern is exactly 
reproduced in every neighbouring cells. Because the loops are staggered within each unit cell (fig. 2.B), there is no 
net magnetization contrary to a ferromagnetic phase. Therefore, the magnetic state corresponds to a Q=0 orbital 
AF order (Q=:0, AFO). 

Motivated by this approach, we developed a neutron diffraction experiment to evidence such a magnetic order. 
Fig. 2.C represents the magnetic intensity in the reciprocal space of a quadratic Cu02 plane taking the magnetic 
structure given in Fig. 2.B. As the state preserves the TSL, the magnetic scattering is only expected on top of 
nuclear Bragg peaks for integer values of H and K. Importantly, due to the existence of two staggered moments 
in each unit cell, no intensity is predicted for H=K=0. In CC — Bn state, 4 classical domains can exist [11]. 
One CC state is shown in Fig. 2. A. The three other states can be obtained by reversing the current flow or by 
rotating Fig. 2. A by 90^ [32]. Two domains should produce the neutron scattering pattern given in Fig. 2.C and 
the two others should produce a similar scattering pattern but rotated at 90°. The scattering pattern in Fig. 2.D is 
obtained by assuming the same population distribution for each domain, restoring the quadratic C4 symmetry in 
the total scattering pattern. One sees that the most intense magnetic peak would be for (H,K)=(zbl, 0) = (0, ±1). 

2.2. Polarized neutron technique 

To separate the nuclear and magnetic scatterings - superimposed at the same Bragg position - requires polarized 
neutron diffraction experiments. As shown in Fig. 2.E, the nuclear scattering interaction conserves the neutron 
spin whereas the magnetic scattering interaction can flip the neutron spin. In a polarized neutron diffraction 
measurement, the incoming neutrons are in the same spin state and one analyses the spin state of the scattered 
neutrons. On can then distinguish the spin-flip (SF) scattered intensity from the non-spin-fiip (NSF) one. This 
technique is extremely powerful. Unfortunately, the (NSF) nuclear intensity is much larger than the (SF) magnetic 
intensity, typically one thousand time larger. The difficulty of the experiment then resides in the capability of 
producing a high polarized neutron beam, reliable and stable in time and in position. In the case where the magnetic 
order is short range, the magnetic intensity is redistributed in momentum space. In this limit, the technical issue 
is related to the weakness of the magnetic signal/background ratio, but the polarization stability is no more 
relevant. These experimental constraints are the best satisfied using the polarized neutron diffraction on the 4F1 
triple-axis spectrometer around the reactor Orphee at the Laboratoire Leon Brillouin (LLB), Saclay (France) 
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Figure 2. A) Schematic description of the circulation currents flowing inside a Cu02 plaquettes. In the CC — On phase, there are 
two currents loops turning clockwise and anti-clockwise. Reversing the current loops and/or rotating the flgure by 90° allows one 
to build the 4 classical domains of the CC — On phase. B) The current loops generate staggered orbital-like magnetic moments, 
perpendicular to the Cu02 planes. C-D) Neutron scattering patterns for a C) single- and D) a multi-domain CC — On phase, 
respectively. In these two flgures, a typical magnetic form factor has been used to compute the magnetic structure factor. E) 
The CC — On phase is a Q=0 AF phase (which belongs to the Eu irreductible representation of the point group [30,31]): it 
breaks time reversal symmetry, but preserves lattice translation invariance. As a consequence, the nuclear and magnetic intensity 
are superimposed on Bragg reflections. The polarized neutron scattering technique allows one to disentangle nuclear and magnetic 
scatterings, which appear selectively in the non-spin-flip (NSF) or spin-flip (SF) channels. 

[33,34,35,36]. The observation of a magnetic scattering whatever the experimental difficulties is a significant sign 
of the universality of the phenomenon. 

As discussed in [33,34,35,36], the polarized neutron diffraction setup is similar to that originally described in 
[37] with longitudinal polarization analysis (see also refs. [38,39] in the context of high-Tc cuprates). A polarized 
incident neutron beam is obtained with a polarizing mirror (bender), the neutron energy is fixed at Ei = 13.7 meV. 
The polarization analysis is performed with an Heusler analyser. A radio-frequency Mezei flipper is put between 
the bender and the sample in order to reverse the incident neutron polarisation. A pyrolythic graphite filter has 
been put before the bender to remove higher order neutrons. The direction of the neutron spin polarization, P, 
at the sample position is selected by a small guide field H of the order of 10 G. 

In polarized neutron scattering technique, one can define the inverse of the flipping ratio (R~^) as = 
Isf/Insf where Isf and Insf stand for the SF and NSF neutron intensities, respectively. For a perfectly 
polarized neutron beam, should be equal to in absence of any magnetic scattering. In real experiments, 
this is not the case. As a results of the unperfect neutron beam polarization, there is a leakage of the NSF 
intensity into the SF channel. At the spectrometer 4F1, a highly polarized neutron beam is routinely obtained 
with ranging between 1/40 and 1/100. In absence of a magnetic order, R~^ exhibits a smooth behaviour with 
temperature, ideally constant but in reality one typically observes a slight drift of R~^ upon cooling. The origin of 
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this temperature dependence is unclear and varies from one experiment to another. It is hkely due to imperfections 
in the neutron polarization and the displacement of the sample in the neutron beam when changing temperature. 
For a better data analysis, one needs to care about this change of the flipping ratio with temperature. R{T) is 
measured from the evolution of flipping ratios at Bragg peaks where the magnetic order is not present. In the 
presence of an additional magnetic scattering, an upturn at low temperature in shows up. We prove that 
method to be eflicient enough to see weak magnetic moments on top of nuclear Bragg peaks, see e.g. 

the determination of the A- type antiferromagnetism in Na cobaltate systems [40]. 

In an unpolarized neutron diffraction measurement and for a magnetism associated with unpaired electrons, 
the interaction between the neutron spin and the magnetic moments in the sample is of dipolar type [41]. As 
a consequence, only the components of the order moment M which are perpendicular to Q the transferred 
momentum contribute to the magnetic scattering |Fm|^ [41]: \Fm\^ oc |M_l|^ with M_L = M - (M.Q)Q with Q 
the unit vector along Q (|Q| = 1). In a polarized neutron diffraction, the amount of magnetic scattering in the NSF 
channel is proportional to |M^.P|^, the neutron polarization vector being normalized (|P|^ = 1). The rest of the 
magnetic scattering (i.e |M^|^ — |M^.P|^) appears in the SF channel. Thus, for the neutron spin polarization 
P//Q, the full magnetic scattering appears in the SF channel. 

2.3. Samples 

Experiments were performed on about a dozen samples of either the YBCO system [33,34], or the Hgl201 
system [35]. These systems, whose crystal structures are given in Fig. 4, allow us to cover a large range of doping, 
especially in the underdoped regime where the pseudogap occurs. In both systems, it is worth pointing out that 
the experimental results are very consistent with each other and with the expected phase diagram, whatever the 
origin of samples, excluding extrinsic effects caused by impurity phases as it has recently argued from a single 
sample [42]. In Hgl201 system, this study has been made possible due to the recent breakthrough in single crystal 
synthesis of samples with mass 1 g [43]. 

The scattering Bragg wave vector Q=(H,K,L) is given in units of the reciprocal lattice vectors, a* 6* = 27r/a 
and c* = 27r/c. Most of the data have been obtained in a scattering plane where all Bragg peaks like Q=(H,0,L) 
were accessible. In the YBCO, the crystal structure is orthorhombic and plain samples are twinned such that 
(H,0,L) = (0,K,L). However, by applying an uniaxial pressure at high temperature (see e.g. [19]), it is possible to 
detwin YBCO single crystals. In such a case, one can distinguish scattering along a* and b*. 



3. Magnetic order in the pseudogap phase 

3.1. Magnetic order in YBCO and Hgl201 

As discussed above, the magnetic scattering is expected to be the largest at in-plane Bragg indices (H,K) = 
(±1,0) = (0,±1) for any integer L value along c*. To evidence the magnetic signal, we need to look for Bragg 
peaks where the structural nuclear intensity is reduced. In all cuprates, this situation is better realized for L=l. The 
Bragg peak intensities in both NSF and SF channels are shown in Fig. 3. A for a detwinned YBCOe.e sample[33]. 
The NSF intensity slightly increases with decreasing temperature as it is expected at small |Q| due to the Debye- 
Waller factor. In contrast, the SF intensity departs from this behaviour below Tmag=220 K evidencing a magnetic 
scattering \Fm\^ developing on top of a smooth background in the SF channel, given by the polarization leakage 
from the nuclear intensity \Fn\^ /R- From both intensities, one can estimate a ratio \Fn\^ /\Fm\^ of ^400 for this 
detwinned YBCOe.e sample. Typically, this ratio is larger in twinned YBCO and Hgl201 samples. Converting the 
relations given in Fig. 3, one can deduce the magnetic scattering cross-section as it has been performed for many 
samples (Fig. 3.B) corresponding to different doping levels ranging from an underdoped YBCOe.s sample to an 
overdoped (Y,Ca)BCO sample [33]. These data were calibrated in absolute units by scaling the Bragg intensity 
to the strong nuclear Bragg peak Q= (0,0,4). One sees in Fig. 3.B that the magnetic scattering appears below 
a temperature, Tmag- Both Tmag and the magnitude of the magnetic scattering increase with decreasing hole 
doping as expected for the pseudogap phenomenon. At variance, similar analysis at Q= (0,0,2) (empty symbols 
on Fig. 3.B) show no additional magnetic scattering below Tmag in agreement with the expected null magnetic 
structure factor for orbital moments (Fig. 2.D) related to a CC order[9,10,ll]. At large |Q|, the magnetic form 
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Figure 3. A) YBCOe.e^ Temperature dependence of the SF intensity (red) and the NSF intensity (blue) measured at Q=(0,1,1) for 
a polarization P// Q. The formula give the expected scattering in both channels where l-Fivl^ and |FmP are the nuclear and the 
magnetic scattering respectively. Both curves in A were rescaled at high temperature by dividing the NSF intensity by R=40 (the 
flipping ratio for that experiment). B) (Y, Ca)BC06+x: Temperature dependence of the magnetic intensity for P//Q at Q=(0,1,1) 
(= (1,0,1)) (full symbols) and Q=(0,0,2) (open symbols). From top to bottom, the hole doping increases through the gradual change 
of the oxygen stoichiometry or extra substitution of Ca for Y. (adapted from [33]). 



factor is expected to considerably reduce the signal. Accordingly, measurements at the Q= (2,0,1) reflection show 
no magnetic scattering [34]. 



3.1.1. Magnetic moment 

Next, from |Fm|^ (Fig- 3.B), one can estimate the magnitude of the ordered magnetic moment |M| defined using 
the magnetic neutron cross section for the CC — Gn phase. Keeping in mind that there are 2 opposite orbital 
moments per Cu02 plaquette, the neutron cross section of the magnetic order given in Fig. 2.B and for a Bragg 
position Q=(H,0,L) simply reads 

|Fm|' = rlf{Qfp{Lf4sm\27rxoH)\M^f (1) 

where ro=-0.54 10""^^ cm {i.e. Tq = 290 mbarns) is the neutron magnetic scattering length, xo = 0.146 is the 
position of the magnetic moment within the unit cell {i.e. the center of the triangle in Fig. 2. A). 

f{Q) stands for the magnetic form factor comprised by principle between and 1. The calculation of the 
magnetic form factor necessitates in principle a detailed knowledge of the real space extension of the orbital 
moments, involving the Cu-da,2_^2 and 0-px,y orbitals. So far, such a calculation has not been performed for the 
observed magnetic order. We take an arbitrary estimate of f{Q)^ — 0.5 around Q=(1,0,1). This crude assumption 
has been previously made for YBCO and Hgl201 [33,34,35] to be able to make a first estimate of the magnetic 
moment. 

P{L) is the magnetic structure factor along c* related to the magnetic structure within the (Cu02)2 bilayer (in 
the case of the YBCO system). As shown from the L-dependence on three different peaks [28,33,34], the coupling 
along c* between the magnetic moments has to be ferromagnetic as the largest magnetic intensity is observed for 
L=0. Accordingly, a natural form for P{L) is then 2cos(7rzL) with 2;=0.29 is the reduced distance between Cu02 
layers in YBCO [28]. 

Considering, on the one hand, a magnetic moment M with the 3 components {Ma, Mb, Mc) along the main 
crystallographic directions and, on the other hand, the existence of 4 domains in the CC — On phase (built from 
Fig. 2. A), the square of the magnetic components perpendicular to Q reads: 
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Hole doping p 

Figure 4. A) Crystal structure of YBa2Cu306+x (a=3.85 A, c=11.7 A). It contains 2 Cu02 planes per unit cell and is orthorhombic. 
Cu sites are embedded into oxygen pyramids and there is a dimpling of Cu02 planes. B) Crystal structure of HgBa2Cu04+5 (a=3.87 
A, c=9.5 A). There is only one Cu02 plane per unit cell and the system is tetragonal. Cu sites are located inside oxygen octahedra. 
Copper and planar oxygens lie exactly in the same plane. C) Hole doping dependence of the superconducting critical temperature 
Tc (stars) and the ordering temperature T^ag of the (Q=0, AFO) magnetic state (bullet). Red and light blue symbols correspond 
to YBa2Cu306+x (YBCO) and HgBa2Cu04-(-5 (Hgl201) systems. The relationship between p and Tc has been established by 
a systematic study of the c lattice parameter in YBCO where the true doping can be simply estimated in the oxygen-ordered 
orthorhombic phases[44]. Unfortunately, the exact doping level in Hgl201 cannot be determined with accuracy. However, the Tc 
curve exhibits the same shape in doping as YBCO giving confidence that the same relationship between p and Tc can be used. The 
size of symbols is proportional to the magnitude of the ordered magnetic moment, (adapted from [35]). 

iMxl^ = ^[1 + C-^fWl, + [1 - C-^fWl (2) 

with Ma,h — \/ + the magnitude of the in-plane magnetic component. At Q=( 1,0,1), {^^Y is equal to 
0.1. If M were perpendicular to the Cu02 plane, then |M^|^ would be '-^ |M|^. As will be shown in the next 
section, M is not simply along the c axis, but rather tilted at ^ 45 ^. One gets a rough estimate |M^|^ ^ ||M|^. 

Finally, using these various approximations, one obtains a magnetic moment of ^ 0.1/is when \Fm\^ ^ 1 mbarn 
in Fig. 3.B. This value decreases when increasing hole doping. It should be emphasized that this moment is given 
per triangular current loop in the CC — Gn model of Fig. 2. A and corresponds to the order of magnitude expected 
theoretically [9]. Of course, the precise value M would change depending on both the model and the accuracy of 
the different parameters appearing in Eq. 1. However, the order of magnitude of the ordered moment is actually 
quite robust whenever other assumptions are made. In Hgl201[35], the same expression Eq. 1 can be used with 
/3(L) = 1 as there is a single layer per unit cell. Using similar assumptions, the magnetic moment per triangle can 
reach ^0.2 iib in the strongly underdoped state[35]. 

3.1.2. A true phase transition 

We then identify this magnetic order to the hidden order parameter of the pseudogap phase in YBCO [33, 34]. 
Similar measurements in Hgl201[35] display the same doping dependence which can be summarized in Fig. 4.C 
where the doping level is determined for both YBCO and Hgl201 using the same updated relationship between 
the doping and Tc [44] . It should be emphasized that the temperature Tmag matches the pseudogap temperature 
T* deduced from resistivity measurements [33, 35]. For both systems, the magnetic ordering temperature Tmag 
extrapolates to zero at p=0.19, which matches the expected end point of the pseudogap phase as deduced from 
various physical properties such as magnetic susceptibility, entropy or resistivity [5]. This points towards the 
existence of a quantum critical point at p=0.19. Clearly, the neutron data provides a strong support in favor 
of a true phase transition at T* [3,45]. In addition, high resolution magneto-optic (Kerr effect) measurements 
consolidate this conclusion as it also evidences a time reversal breaking symmetry in the pseudogap phase of 
YBCO [46]. 

The pseudogap affects all physical properties. However, the value of T* for a given doping varies through the 
literature depending on the data analysis or on each experimental technique. That is actually the reason why T* 
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Figure 5. A) Schematic description of the three orthogonal neutron polarizations (see text). The shaded area represents the scattering 
plane defined by the directions (100) and (001). B-E) Temperature dependencies of the intensities in the SF channel (magenta): B) 
for P // Q at Q =(1,0,1) and (C) at Q=(2,0,l), (D) for P // z and (E) P ± Q at Q=(1,0,1). For all panels, blue symbols indicate 
the polarization leakage, given by the NSF intensity divided by the temperature dependence of the fiipping ratio, R(T). R(T) is 
found by a fit in panel C of the ratio NSF/SF of the Bragg peak Q=(2,0,l) as R(T)=R(300 K){l+0.02(l-T/300) F) Orientation of 
the magnetic moment with respect to the c axis (green arrow), as deduced from the polarization analysis (adapted from [34]). 

has been largely considered as a crossover phenomenon. In order to compare the effect of the pseudogap on a given 
physical properties, one needs to know how this quantity can couple to the order parameter of the pseudogap 
phase. For instance, it has been argued [47] that the CC — Gn theory could be mapped onto the Ashkin- Teller 
model (a model with a pair of Ising spins at each site which interact with neighboring spins through pair-wise and 
four-spin interactions). Monte Carlo simulations [45] of the Ashkin- Teller model show no sharp thermodynamic 
anomaly at the phase transition for a large set of parameters. That may corresponds to that observed in cuprates 
at T*[5]. Furthermore, the order parameter exponent /3 =0.18 [34] that one can extract from the neutron data 
falls into the proper range for such an Ising type model. Another example of a thermodynamic signature of the 
existence of a phase transition in the pseudogap state of underdoped YBCO is given by recent high-precision 
magnetization measurements [48]. The temperature derivative of the uniform susceptibility indicates a singular 
point at a temperature corresponding to Tmag- This can be understood [45,48] as a biquadratic coupling of the 
magnetic order parameter with the uniform magnetization in a way similar to what happens in antiferromagnet. 



3.1.3. Polarization analysis 

Polarized neutron diffraction technique allows us to determine the direction of the magnetic moment through 
a polarization analysis[37]. As the interaction of the neutron spin with magnetic moments is of dipolar type, only 
the magnetic components perpendicular to the wavevector Q are measurable [41]. Furthermore, only a magnetic 
component perpendicular to the neutron spin polarization P is spin-flip whereas a component along P is non- 
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Figure 6. YBCOe.e: A) [100]/[001] scattering plane and schematic description of a scan along the [001] direction across Q=(1,0,1). 
B) Difference between L-scans at 70 K and 300 K in the SF channel and for polarisation P//Q (magenta). The magnetic signal, 
centered at Q=(1,0,1), is resolution limited. The blue symbols stand for the resolution limited nuclear Bragg scattering at Q 
= (1,0,1) in the NSF channel (adapted from [34]). 



spin- flip [41]. The amount of magnetic scattering in the SF channel therefore depends on the choice of P. 

For a polarized neutron analysis, one usually introduces a set of three orthogonal directions: The x axis is 
parallel to Q, the y axis is perpendicular to Q in the scattering plane and the z axis is perpendicular to the 
scattering plane. For a neutron spin polarization along x (P//Q = H^), the full magnetic intensity appears in 
the SF channel: /p//Q oc M]_y + M^^, M±y and M±z are the components of along y and z respectively. For 
a neutron spin polarization applied along y (P-LQ = H^), only one part of the magnetic intensity appears in the 
SF channel: Ip±q oc M]_^. The complementary magnetic intensity /p//z ^ ^±y shows up in the SF channel when 
the neutron spin polarization is applied along z (P//z = H^). As a result, in the SF channel, one can deduce 
a specific polarization sum rule for the neutron intensity: /p//Q = Ip±Q + ^p//z- It is worth noticing that this 
polarization sum rule is only fulfilled for a magnetic scattering in the absence of chirality in the system. 

The polarization analysis has been performed in YBCO[33,34] for the Bragg peak Q=(1,0,1) in twinned samples 
as shown in Fig. 5. In first approximation, the flipping ratio is constant in temperature [33, 35]. To improve the 
data analysis, one needs further to fit it empirically by a straight line [34,49] as it has been done in Fig. 5.C. 
First, the magnetic scattering is larger for P//Q as it should be. Next, the polarization sum-rule above is fulfilled 
proving the magnetic nature of the signal observed below Tmag- The same conclusion has been also demonstrated 
in another YBCOe.e sample [33] as well in recent data obtained in Hgl201 [49]. 

The data reported in Fig. 5 were measured in the scattering plane given by the directions (100) and (001): 
the z axis therefore corresponds to the (010) direction. The orbital magnetic moment in a CC-phase is then 
expected to be parallel to c (Fig. 2.B). Thus, one should expect /p//Q = ^p//z and /p±Q=0. At variance, one 
observes a magnetic contribution for all neutron polarizations (Fig. 5.B,D,E). More specifically, /p//Q=1.4 mbarn 
and Ip±Q = /p//z=0.7 mbarn in the YBCOe.e sample[34] shown in Fig. 5. Qualitatively, the fact we observe a 
non-zero magnetic intensity for P//z (Fig. 5.D) proves that a large part of the magnetic moment is pointing along 
c. However, the intensity for P//z is smaller than the one for P//Q indicating that an in-plane component is also 
present. 

Next, to discuss the magnetic moment direction, we use the same simple model used in previous sections to 
evaluate the order of magnitude of the ordered magnetic moment. We assume a model with collinear moments 
pointing along a generic direction given by, M=(Ma, M^, Mc) with no preferential in-plane direction (Fig. 5.F). One 
then needs to calculate the neutron intensity from Eq. 1 and Eq. 2 for Q=(1,0,1) and for each neutron polarization 
[41]. That gives Ip±q oc [|Mf 5], /p//^ oc [0.05Mf ;, + 0.9Mc ] and the sum of both terms for P//Q. One can define 

the angle, between the magnetic moment and c (see Fig. 5.F) which can be written as tan((/)) = ^/ ^/M^ . 
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Figure 7. A) Body centered crystal structure of monolayer system La2-xSrxCu04. The hole doping is achieved by substitution 
of Sr^"^ for La"^"*" and the hole doping (p) is equal to the Sr content (x). In the low temperature orthorhombic (LTO) phase, the 
orthorhombicity appears along the [110] direction in tetragonal notation. B) Phase diagram of La2-xSrxCu04 system as a function 
of hole doping. In the low doping regime below p=0.05, one observes 2 states: the antiferromagnetic (AF) state and the cluster 
glass state with diagonal stripes. In the metallic regime above p=0.05, the system becomes superconducting below Tc and below 
the freezing temperature Tg a cluster glass state with longitudinal stripes can be observed. This phase becomes more stable upon 
approaching the critical hole doping p=l/8 and always tends to develop at the expense of the superconducting state (from [56]) C) 
Cartoon picture showing a bond centered longitudinal stripe state close to p=l/8. 

Using the result of Fig. 5 at Q=(1,0,1), one gets = 55 =b T'^ [34]. At Q=(1,0,0), the angle is found to be smaller 
= 35=1=7'^ [34]. It is unclear at the moment what is the meaning of such a difference. The given errors are related 
to statistical error, however, due to possible systematic bias in the data analysis, the actual error might be larger. 
One can then give a conservative estimate of = 45 di 20^^ valid over few samples and Bragg spots [33] . The recent 
data in Hgl201 leads to the same tilt of the moment relative to the c axis [49]. One then sees that this effect does 
not depend on the number of Cu02 per unit cell but is intrinsic to that novel (Q=0, AFO) magnetic order. The 
implication of this finding will be discussed in the last section. 

To estimate the angle 0, we here assumed a collinear magnetic structure for simplicity. However, it should be 
noticed that more complex magnetic structures might exist, in particular those involving non-collinear moments. In 
such cases, different magnetic diffraction patterns could be obtained for each magnetic component in the neutron 
scattering cross-section. In principle, this might yield a different estimate of the tilted angle. 

3.1.4. Long range order 

So far, only the intensity at the Bragg position has been discussed. One crucial question about this (Q=0, AFO) 
broken symmetry is whether the order is at long range or not. To answer that question, one needs to perform 
scans in the reciprocal space without loosing the polarization condition. At the level of accuracy necessary to seize 
the magnetic peak, this is actually only working properly along c* (sketched by the blue arrow in Fig. 6. A) which 
almost corresponds to a rocking scan of the sample. Fig. 6.B depicts the magnetic scattering, obtained by taking 
a temperature difference of the SF scattering for the P//Q. The blue points represent the non-magnetic NSF 
Bragg scattering scaled to the same amplitude as the magnetic scattering. The q-width of the NSF Bragg peak is 
given by the resolution of the spectrometer: here, the full width at half maximum is 0.013 A-^ The SF and NSF 
intensities superpose very well. This implies that the q-width of magnetic peak is also limited by the spectrometer 
resolution, suggesting that the magnetic state is ordered at long range along c [33,34]. From the resolution width, 
one can determine a lower bound for the magnetic correlation length of 75 A. 

3.2. The case o/ La2-xSrxCu04 

The monolayer system La2-xSrxCu04 (LSCO) is the archetypal high-temperature copper oxides superconduc- 
tors, mostly because by changing continuously the Sr content one can cover all the regimes of the cuprates phase 
diagram from the insulating AF state upto the overdoped state with Fermi-liquid properties. However, the maxi- 
mum SC temperature reaches at most 37 K, a much smaller temperature than in other monolayer system such as 
Hgl201 (97 K). In comparison with systems like YBCO or Hgl201, LSCO displays many drawbacks. As a function 
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of temperature, the system undergoes a structural transition from a high temperature tetragonal phase (HTT) to 
a low temperature orthorhombic phase (LTO). The hole doping is achieved by cationic La^^/Sr^^ substitution. 
Such a substitution reduces the tilt of the oxygen octahedra which is at the origin of the structural distortion. 
The impact of this structural transition in transport or thermodynamic measurements often blurs an accurate 
determination of properties like the pseudogap temperature. Furthermore, the Sr^^ ions, located close to Cu02 
planes, can introduce random electrostatic potentials within the planes. LSCO is therefore a rather disordered 
system: this can be seen in the resistivity measurements, characterized by a large residual resistivity [50] and in 
NMR data, where the NMR linewidth are much larger than in cuprate whose Tc is about 90 K at optimal doping 
[51]. In many cuprates families, an upturn of the resistivity can be observed at low hole doping on cooling down 
[52] . Such an upturn is characterized by a change of sign of the first derivative of resistivity as a function of 
temperature (-^ < 0), as observed in semiconductors or insulators. In addition, when applying a large external 
magnetic field in order to reduce or remove superconductivity, this insulating-like behavior at low temperature 
can be detected up to a critical hole doping pm, above which a metallic behavior is recovered [53]. pm is not 
a universal critical hole doping among cuprates: in LSCO, pm is as large as 0.16, corresponding to the optimal 
doping. In less disordered materials, such as YBCO, pm remains smaller than 0.09 [24]. The peculiarities of LSCO 
system appear also in characterization of the pseudogap state. While in YBCO, T* can be easily determined 
from resistivity measurements or NMR knight shift, the characteristics of the pseudogap in LSCO are typically 
ill-defined [4]. 

Beyond, its intrinsic disorder, LSCO system is characterized by a pretty strong tendency towards the formation 
of inhomogeneous spin and charge distributions in the Cu02 planes, usually associated with ID stripes [13]. Stripes 
can be viewed as filamentary charge organisations, separating AF domains in anti-phase. A picture of bond centered 
stripes is given in Fig. 7.C for a doping of 1/8, doped holes are confined into 2-leg ladders, separating hole poor 
domains with an antiferromagnetic arrangement of S=l/2 Cu spins. The 2-leg ladders with doped holes behave 
as antiphase boundaries between hole poor AF spin stripes. In LSCO, the low energy spin excitation spectrum is 
dominated by strong incommensurate (IC) spin fiuctuations around the planar antiferromagnetic (AF) wavevector 
at Qic = Qaf =b ((^,0) = Qaf =b (0,5) [54] with Qaf = (0.5,0.5). These low energy IC fiuctuations are usually 
interpreted as the fingerprint of fiuctuating stripes. Local probes (NMR, /iSR,) [55,56] have established that these 
fiuctuations can become static at very low temperature. These observations suggest that fiuctuating stripes can 
either freeze at low doping or can also be pinned down by disorder, yielding a glassy static stripe phase. As shown 
in Fig. 7.B, this tendency towards stabilization of static stripes is reinforced around a doping of 1/8 where stripes 
are better pinned to the lattice. In isostructural compounds La2-xBaxCu04 [57] and (La, Nd)2-x(Sr, Ba)xCu04 
[18], where superconductivity is strongly reduced, the situation is even clearer as a spin density wave (SDW) order 
develops at Qic around 30 K, in addition to a charge density wave (CDW) at 2Qic at higher temperature [18]. 

LSCO is then a very interesting and complex system at the crossroad of many different electronic instabilities. 
What might be the Q=0 AFO phase in that system has to be answered. Therefore, we looked recently for the Q=0 
magnetic order in a LSCO sample with x=0.085 with Tc=22 K [36]. 



3.2.1. Absence of long range order in LSCO 

Similar measurements as YBCO and Hgl201 have been repeated in LSCO. The magnetic moments of the Q=0 
AFO always scatter neutrons at the same positions as the Bragg peaks of the crystallographic structure, but in 
contrast to YBCO and Hgl201, the crystal structure of LSCO is not primitive with a stacking along c with CuOe 
octahedra shifted along the lattice diagonal (Fig. 7. A). Resulting from the body centered tetragonal phase at high 
temperature, the low temperature structural phase is face centered orthorhombic: Bragg positions like (1,0, L) with 
even L are forbidden. The search of the long range Q=0 magnetic has therefore been carried at the Bragg point 
Q=(1,0,1). As in other cuprates [33,34,35], this Bragg position offers the best compromise for a weaker nuclear 
Bragg peak having the proper symmetry for the Q=0 AFO. A long range 3D Q=0 magnetic order would give rise 
to an upturn of the inverse fiipping ratio {R~^) on cooling down. In Fig. 8, we report on the same scale for 
LSCO and YBCOe.e [33] for which a magnetic signal is clearly visible below 220 K. The temperature dependence 
of for LSCO only shows some fiuctuations due to minute changes of the polarization of the spectrometer. 
Therefore, the orbital-like long range order present in other cuprates is either absent in LSCO or too weak to be 
experimentally detected. As nuclear Bragg intensities in LSCO and YBCO are similar and using the error bars 
of Fig. 8, one can give an upper estimate of less than 0.02 /is for a 3D ordered orbital-like moment in LSCO, 
instead of ~ 0.1 /is in YBCOe.e for a hole doping p~ 0.12 [33]. 
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Figure 8. Temperature dependencies of the inverse of the flipping ration (-R~^) measured at Q=(1,0,1) for P // Q : for YBa2Cu306.6 
(YBaCuO - p=0.12 - green symbols from [33]) and La2-xSrxCu04 LaSrCuO - x=0.085 - red symbols). The upturn of below 
^mag =220 K indicates the appearance of the 3D Q=0 AF order upon cooling down. This feature is absent in T-dependence for 
La2-xSrxCu04 system, pointing out the absence of the 3D Q=0 AF order. 

3.2.2. Short range and 2D orbital-like order 

However, there are other possibiUties for Q=0 AFO as shown in Fig. 9.C: The magnetic scattering can be 
either present at forbidden Bragg points for the structure, or even exist at any position along c*. We checked 
experimentally these two scenarios [36]. We first studied the magnetic scattering around the forbidden Bragg 
reflection (1,0,0). Fig. 9. A shows a scan in the spin-flip channel for P//Q along the diagonal [Irl^O] direction 
at 30 K and 120 K. A broad peak is observed there although no peak is seen in the NSF channel [36]. The 
observation of such a magnetic intensity at a forbidden Bragg peak may mean that the stacking of the orbital- 
like moments does not preserve the lattice symmetry. The low temperature scan can be fitted by a Gaussian 
peak on top of a sloppy background. From the deduced peak width Aq = 0.11 ± 0.02 r.l.u , one obtains a short 
correlation length of ^no = 2/Ag 11 ± 2A after resolution deconvolution. The magnetic scattering does not 
appear at the forbidden Bragg positions only. Indeed, additional scans at non-integer L exhibit a similar broad 
peak independently of L [36]. For instance, we report on Fig. 9.B the scan for L=0.45 along the [1,0,0] direction 
in which the signal appears to be broader along that direction yielding ^loo = 8 i 3 A. It should emphasized that 
a similar signal is also expected at the Bragg position L=l but the signal is too small to be detected even with 
polarized neutron diffraction. In contrast to YBCO and Hgl201, the magnetic signal in LSCO is then at short 
range and bidimensional (2D) as sketched by the yellow area in Fig. 9.C. but it is still localized around the same 
planar wavevectors respecting the translation symmetry of the Cu02 plane, say Q2i:) = (l,0), corresponding to the 
Q=0 AFO. 

Fig. 9.D depicts the polarization analysis of the short range magnetic order in LSCO. As in YBCO and Hgl201, 
the magnetic scattering occurs for all neutron polarization directions but is weaker than for P//Q, fulfilling the 
above-mentioned sum rule. The fact that /p±q ^ ^p//z means that the moments are pointing in any preferential 
direction. This can be explained by considering either a rather disordered magnetic state or a similar situation as 
in the case of YBCO [33,34] and Hgl201 [35,49], where the moments are tilted from the c* axis by ~ 45^^. 

The reported 2D magnetic intensity is rather weak with a ratio signal/background of about only 5% [36]. One 
then can ask whether this magnetic intensity is sufficient to produce a large enough local moment. For a 3D long 
range magnetic order, the magnetic intensity at a magnetic Bragg position is directly proportional to square of the 
magnetic moment (Eq. 1). For a 2D short magnetic signal, the magnetic intensity is redistributed in momentum 
space. In that case, this is the integrated magnetic intensity over the momentum space that scales with the square 
of the local magnetic moment. At first, one needs to put the 2D magnetic intensity, Imag, in absolute units after 
taking into account the spectrometer resolution. Imag is described as 



ImagiQ) — Iq exp 
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Figure 9. Evidence for a short range 2D Q=0 AF order in La2-xSrxCu04: A) Diagonal scans across Q=(1,0,0) in the SF channel 
for P//Q. B) Longitudinal scans across Q=(l,0,0.45) in the SF channel for P//Q. C) Schematic description of a longitudinal scan 
across the magnetic rod (1,0, L) in the [100]/[001] scattering plane. At Low temperature, scans show the appearance of a signal 
centered at Q=(1,0,L) on top of a smooth background. The signal exhibits a Gaussian profile, whose intrinsic full width at half 
maximum (Ag) allows an estimation of the magnetic correlation length ^ (=2/Aq). D) Polarized neutron scattering analysis showing 
the magnetic intensity in the SF channel for P//Q (Hpx - red squares), P±Q (Hpy - green bullets) and P//z (Hpz - blue triangles). 
The black symbols stand for the SF background (BG), which does not depend on the neutron spin polarization. E/ Temperature 
dependence of the magnetic intensity, calibrated in absolute units (mbarns) (adapted from [36]). 

where ^ is the magnetic correlation length assumed to be isotropic within the (H,K) plane for simplicity. Q2z:>=(l,0) 
is the planar wavevector around which the magnetic scattering occurs. Imag is L-independent and assumed to be 
static at the neutron energy scale, meaning with a characteristic energy lower than ^ 1 meV. The 2D magnetic 
intensity, Imag, is then convolved by the Gaussian resolution of the spectrometer. Its amplitude is further normal- 
ized in absolute units by comparison with the nuclear Bragg peak Q=( 1,0,1). Then the Q-integrated magnetic 
intensity, Smag is obtained by integration of Eq. 3 as Smag = / d^QImag{Q) / f d^Q — (a/^)^/o/(47r ln(2)) where 
a=3.82 A is the planar lattice parameter. Taking a mean value of ^=10 A and assuming that it is temperature 
independent, one obtains the integrated magnetic intensity, Smag- As shown in Fig. 9.E, Smag saturates at 1.2 
mbarns at the lowest temperature, a value equivalent to the Bragg magnetic peak intensity in YBCO (Fig. 3.B) 
[33,34]. 

Next, one can determine the local magnetic moment Mioc following the magnetic neutron cross section (Eq. 1) 
applied to the case of LSCO meaning i) |M^|^ = 2/?>Mf^^ corresponding to a distribution of static magnetic 
clusters without any preferential orientation of magnetic moments from one magnetic domain to another , and 
ii) I3{L) — 1 because LSCO has a single Cu02 per unit cell. Using the same mere assumptions as for YBCO and 
Hgl201, the magnetic local moment Mioc ^ OA/j^b is obtained, remarkably, of similar amplitude as the long range 
order Q=0 AFO in both YBCO [33,34] and Hgl201 [35]. 

The temperature dependence of this integrated magnetic intensity, displayed in Fig. 9.E [36], shows a transition 
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Tmag around 120K, in agreement with the absence of a signal in high temperature scans. In YBCO and Hgl201, 
Tmag is associated with the pseudogap temperature T* [33,35]. According to the usual phase diagram of cuprates 
(Fig. 1), a much larger value for T* is in principle expected. However, one the one hand, the pseudogap properties 
in LSCO is less accurate than in other cuprates[4]: that could be related to the observed 2D state, characterized 
by finite correlation length. On the other hand, the sample doping is quite low (8.5 %): the determination of the 
pseudogap behaviour is hampered by the proximity of the metal-insulating crossover as well as by the existence 
of strong AF fluctuations [54]. Actually, several anomalies have been reported close to Tmag for LSCO samples in 
the same doping range: these anomalies appear in the specific heat [58,59], the uniform spin susceptibility [58] and 
the Nernst effect [60] . However, the exact relation between these various anomalies with the short range magnetic 
order is not settled yet. 



4. Discussion 

In three different cuprates families, we observe a novel magnetic order in the pseudogap state. Its order parameter 
breaks time reversal symmetry in agreement with the observation of spontaneous Kerr effect in YBCO [46] and 
of dichroism in the ARPES spectra in another cuprate, Bi2212 [29]. Likewise this novel magnetic order preserves 
the translation symmetry of the lattice, being then categorized as a Q=0 antiferromagnetic order. What are the 
implications of this symmetry ? Independently of any theory, it proves that at least two antiparallel magnetic 
moments per unit cell are necessary to model our observations. Cu spins, which are usually invoked to explain 
magnetism in cuprates, cannot account for the novel magnetic phase as there is only Cu spin per unit cell. 
Its magnetic ordering gives rise to additional Bragg peaks (typically incommensurate) near the (1/2,1/2) wave 
vector[18,19,20]. Our finding implies that one needs to reconsider the role of oxygen orbitals, especially the in-plane 
p-orbitals. Theoretically, that questions the standard effective one-band Hubbard model generally considered to 
encompass the physics of high- Tc cuprates. A priori, our observation implies that a proper starting point has to 
be at minimum a three-bands Hubbard model. 

4.1. Orbital or spin moments ? 

The observation of a Q=0 order in the pseudogap state of cuprates has fundamental implications. In order 
to account of such novel magnetic state, a decoration of the unit cell is necessary with at least two opposite 
magnetic moments within each squared Cu02 plaquette. Indeed, together with the absence of magnetic moments 
in magnetization measurements [48], the zero intensity for H=K=0 ruled out ferromagnetism as the origin of 
the magnetic structure. Our studies therefore suggest that a single Cu02 plaquette contains internal degrees of 
freedom, which have been ignored in most of theoretical models developed for cuprates. This viewpoint has been 
recently experimentally comforted by the observation of electronic nematicity in Bi2212 by STM [61] which also 
suggests intra unit cell electronic patterns. 

Of course, the observation of the novel Q=0 magnetic state rises many questions. What can be the origin of this 
magnetism which preserves the translation symmetry of the lattice? Are we dealing orbital-like or spin moments 
? Unfortunately, the information brought by the neutron experiment is rather scarce in terms of discriminating 
among models as, so far, the magnetic intensity has been only observed at three different Bragg peaks Q=(1,0,L) 
(L=0,1 and 2) and is absent at Q=(0,0,2) and Q=(2,0,l). It is therefore delicate to determine a magnetic structure 
based on that. 

A way to distinguish among orbital and spin moments is to measure the magnetic form factor (/(Q) in Eq. 
1). Indeed, orbital-like magnetic moments, generated by current loops, are spatially extended objects and f{Q) 
is expected to drop down at large |Q| much fast that f{Q) for spins. In the latter case, f{Q) is given by Fourier 
transform of the electronic distribution around magnetic ions. Considering the small number of studied Bragg 
reflections, the interesting approach has turned to a dead-end. For instance, orbital and spin magnetic form factors 
would give no intensity or an extremely weak intensity at Q= (2,0,1), as observed experimentally. To date, there 
is not enough data to tell whether the magnetic form factor is more consistent with orbital moments or spin 
moments. 

Another way to get a deeper insight into the intrinsic nature of the novel magnetic order is to perform an analysis 
of the neutron intensity in terms of the magnetic group symmetry [30, 31]. This analysis should be done relative 
to the crystalline group of cuprates. However, one additional difficulty is here that we even do not know which 
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crystallographic sites can carry the magnetic moments. For instance, in orbital moments model, the moments 
do not sit at given atomic positions. Further, the differentiation of models is as usual complicated due to the 
possibility of multiple domains. All of that implies that the symmetry analysis does not allow to discriminate 
different models, as the ones shown in Fig. 10, based only on the limited neutron diffraction data although these 
models are clearly breaking different symmetries. 

One can nevertheless discuss the allowed symmetries as it has been done by Kaur and Agterberg [30] when 
considering the additional observation of dichroism in ARPES in Bi2212 [29] which, a priori, is related to the same 
broken symmetry phase. The observation of dichroism at the M=(7r,0) implies that both the fourfold rotation 
about z axis (C4) and the diagonal mirror planes {ada, cFdh) are no longer symmetry operations in the pseudogap 
state. That reduces the number of possible magnetic point groups of the pseudogap phase to only two [30]. Further, 
the photoemission matrix elements analysis of the dichroic signal [10] reduces the number of magnetic point groups 
to only one referred as mmm, Eu irreductible representation of the D4h point group [30]. The orbit als moments 
model created by current loops in the CC — Gn state (shown in Fig. 10. A and FIG 2), proposed by Chandra 
Varma and coworkers[10,ll,31], belongs to that magnetic group. 

As a matter of fact, that CC — Gn model gives magnetic intensity at the proper Bragg reflections. In contrast, our 
observation of a signal at the in-plane wave vector (1,0) or (0,1) dismisses the CC — Gi phase (of Big irreductible 
representation) with four current loops per unit cell [9,10] as no signal is expected at these Bragg positions for 
that model. The same CC — Gi phase was also ruled out [10] in Bi2212 by quantitative analysis of the ARPES- 
dichroism [29], only compatible with the CC — Gn state. In principle, other magnetic structures could be also 
compatible with the observed magnetic intensity. One interesting example is a staggered orbital current phase 
proposed by Stanescu and Philipps [62] that we have reproduced in Fig. 10. B. It has actually the same symmetry 
as the CC — Gn but does not explicit the role of oxygens atoms. As a result, it is also fully compatible with the 
magnetic neutron intensity observed around Q=(1,0,L) Bragg spots and its absence at Q=(0,0,L). These authors 
claim that it is possible to generate such a translationally invariant current pattern within a one-band model 
[62] once oxygen orbit als are integrated out. The next nearest hopping, t', across the unit cell diagonal would be 
essential to get the pattern. Other staggered orbital current patterns are also compatible with the observed Bragg 
peaks once the apical oxygen orbit als are taken into account: variational Monte- Carlo calculations in the extended 
Hubbard model indeed show that orbital moments develop on the copper oxygen octahedra [63]. Currents are 
flowing on certain faces of the CuOe octahedra giving rise to orbital moments tilted from c (Fig. 10. C). 

It is worth mentioning that spin based magnetic structure can also account for the neutron data. For instance, 
the model shown in Fig. 10. D [33] with spin moments on the in-plane oxygen atoms possesses the necessary 
ingredients: there is not net magnetic within the unit cell and the more intense magnetic Bragg peaks are for 
Q=(1,0,L) = (0,1,L). This order belongs to the Big irreductible representation and preserves spacial inversion 
contrary to the CC — Gn phase. 
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4.2. Magnetic excitations in Hgl201 



Having observed a long range magnetic ordered state, the next step is to search for the magnetic excitations 
associated with that novel magnetic order. By principle, their observation would be an essential issue to validate 
and discriminate theoretical models. We have then performed inelastic polarized and unpolarized neutron exper- 
iments to look for magnetic excitations related to the pseudogap order. Recent measurements in Hgl201 reveal 
such a collective magnetic mode [64] around 54 meV. A peak is observed in the spin-flip channel and absent in the 
non-spin-flip channel in polarized neutron experiments. 

The mode is associated with the Q=0 magnetic order as the mode intensity is going to zero at the temperature, 
Tmag = T* where the magnetic order vanishes and where anomalies occur in resistivity measurements. This has 
been observed for two different samples having very different T* in the underdoped regime. Furthermore the 
magnetic mode exhibits a very weak dispersion over the whole Brillouin zone. This behaviour drastically differs 
from the magnetic fluctuations usually reported in high-Tc cuprates [1] which occur systematically around the 
planar antiferromagnetic (AF) wavevector Qaf=(0.5,0.5). Near the zone center, these AF fluctuations are absent. 
In contrast, the new excitations exist at any wave vector around nearly the same energy. That is what would be 
expected for an order parameter characterized by a discrete symmetry, such as in the case of Ising-like moments. 

This supports the orbital current phase, CC-Gn [11], which can be actually represented by an Ashkin- Teller 
model [65], characterized by two coupled Ising moments. More specifically, three weakly dispersive modes are 
predicted from the CC — Bn model [32]. The two low energy modes are measurable by inelastic neutron scattering 
as they are related to a flip of an orbital moment by AL=1. They can be classically represented by a rotation 
by zb 90^, clockwise and ant i- clockwise, from a loop current domain such as Fig. 10. A. One of these modes 
is likely reported in [64], whereas preliminary measurements indicate another mode at lower energy. The third 
mode corresponds to a flip of both moments in a domain like Fig. 10. A: it then cannot be measured by neutron 
scattering as it is associated with a change of the orbital moment AL=2. The experimental determination of that 
third excitation by other techniques like resonant inelastic X-ray scattering (RIXS) would be a crucial test for the 
theory. More generally, the observation of these weakly dispersive magnetic excitations challenge the theoretical 
description of high-Tc cuprates where magnetism essentially comes from AF correlation between Cu spins. For 
the moment this observation of novel excitations is limited to Hgl201 system, it has to be generalized to other 
cuprates. 

4.3. Direction of moments 

As shown above in section 3.1.3, the observed magnetic structure indicates a moment which does not point 
along a high symmetry direction but instead is tilted from the c by ~ 45^ (Fig. 5.F). It is worth to mention 
that the Q=0 magnetism displays a large component along c, contrary to usual Cu^^ spin magnetism in undoped 
cuprates. That suggests a distinct origin. 

Clearly, the tilted picture is not consistent with the planar circulating currents models (Figs. 10. A and B) where 
the orbital magnetic moments have to be perpendicular to the current loops, i.e. pointing completely along c axis. 
At first sight, it seems more natural for spins to explain a moment pointing in any direction. However, this view 
might be too simple. Even in the CC picture, a few arguments can be given to explain the origin of the tilt. First, 
the proposed phase of circulating currents flowing on the face of CuOe octahedra [63] obviously gives rise to a 
tilted moment with a proper order of magnitude. Second, even within planar orbital currents, corrections to the 
original CC-Bn model can be presented. In YBCO, for instance, the Cu-0 bond is not strictly perpendicular to 
c as copper and oxygen atoms are not in the same plane, yielding a tilt of about 14°. Still in YBCO, an in-plane 
spin component necessarily occurs due to spin-orbit coupling [66] whose amplitude might be measurable with 
neutrons although it would not be a very large contribution. In Hgl201, these two last points are not present due 
to its simple tetragonal lattice symmetry. Nevertheless, even for simple squared lattice, an in-plane component 
can appear in the CC-Bn model due to quantum corrections of the ground state [32]. Indeed, the CC — Bn 
is characterized by 4 distinct states for a given Cu02 plaquette. If each state is considered independently from 
the others, that gives 4 classical domains. Alternatively, once quantum corrections are introduced, the ground 
of the CC — Bii becomes a quantum superposition of these 4 states. In neutron diffraction measurements, what 
is interpreted in terms of a planar component could simply result from an interference phenomenon: in such a 
case, the tilt angle indicates the degree of admixture of the 4 CC states within the ground state. Within this 
approach, the exact amount of the tilt is still under question, as it is directly related to te dispersion of the new 
Q=0 magnetic modes [3 2]. 
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In addition to other symmetries discussed above, it is interesting to note that the existence of the in-plane 
component breaks the Gz mirror plane. In Landau theory of continuous second order phase transition, this indicates 
that each magnetic component is represented by two different order parameters likely having distinct symmetries. 
This means in principle that both order parameters would also exhibit different critical temperatures. 



4.4. Nuclear Magnetic Resonance and muon-spin rotation: silent local probes 

Besides circulating currents, we have mentioned that the Q=0 AFO could also be produced by spins located 
on oxygen sites. There is however a noticeable difference between a spin-based model and an orbital-based model. 
Spins located on oxygen sites should produce a sizeable broadening of the Nuclear Magnetic Resonance (NMR) 
^^O lines which is not observed [67]. At variance, the staggered orbital moment of CC phases [9,10,11,62] would not 
give rise to any measurable signal in NMR experiments. By symmetry, the effects of orbital-like moments cancel 
out for each nucleus site. As a matter of fact, there is no indication in NMR measurements of a static magnetic 
ordering at T*. In underdoped YBCO, there is for instance no broadening of the ^^Y NMR line which actually 
sharpens up at low temperature [51]. Recently, a ^^Y NMR study has been misleadlingly used to point out a lack 
of evidence for orbital-current effects in the cuprate Y2Ba4Cu70i5-5 [68]. This system contains 2 adjacent Cu02 
planes with different hole doping of about ~ 3 % due to different environment. Owing to this charge imbalance, 
the authors of Ref. [68] expected a finite broadening of the ^^Y NMR line for orbital currents phase and they 
did not observe it. Unfortunately, this expectation is not realistic for all the CC phases respecting the lattice 
invariance symmetry discussed here [9,10,11,62] because the magnetic field at the ^^Y site almost cancels out the 
effect as there are a few opposite current loops around each yttrium nucleus. 

In principle, zero-field muon-spin rotation (/xSR) technique is a very powerful to detect very tiny magnetic 
moments. A signature of the Q=0 AFO reported by polarized neutron should be visible in /j,SR measurements. 
Unfortunately, interpretation of /iSR results have been quite contradictory on that topic over the last decade 
[69,70,42,71]. At first, small spontaneous static magnetic fields of electronic origin were reported in YBCO to 
be intimately related to the pseudogap transition [69]. Later on, the same data were interpreted as being due to 
spatial charge inhomogeneities in relation to the CuO chains [70]. Next, a magnetic order was observed but did not 
evolve significantly with hole doping [42]. In only one sample, a second magnetic component was interpreted as due 
to dilute impurity phases [42] . That signal was misleadlingly associated with the magnetic order that we reported 
in neutron diffraction [34] which in contrast to /iSR is observed systematically in all underdoped YBCO samples 
that have been studied [28,33]. The underdoped YBCO sample carrying a dilute impurity phases in the (/iSR) 
study [42] is known to exhibit a magnetic signal around Qaf= (0.5,0.5), easily detected by neutron diffraction 
[72]. It is notorious that such a kind of parasitic AF phase does not exist in all underdoped YBCO samples [73]. 
It then appears natural to ascribe the /xSR impurity phase to the parasitic AF signal contrary to what has been 
claimed in [42]. Another /iSR in LSCO came to the conclusion of the absence of broken time-reversal symmetry 
in the pseudogap state [71]. All these reports then suggest that /iSR probe is silent with regard to the magnetic 
order reported in polarized neutron diffraction as it is also the case for NMR. 

Actually, two independent theoretical calculations [74,75] estimate the screening charge density distribution due 
to a point charge, such as that of a positive muon /x^, placed between the planes of a highly anisotropic layered 
metal. Clearly, a muon is shown to lead to an observable local perturbation. This questions the interpretation of 
muon-spin-rotation experiments in metallic high-temperature superconductors. For the Q=0 AFO, the field at the 
muon site was expected to be of several mT according to the neutron data. This field is actually reduced by more 
than 2 orders of magnitude due to muon related screening [74]. That seems to resolve the apparent contradiction 
between the polarized neutron diffraction and fiSR experiments. 



4.5. Electronic liquid crystal 

It has been recently proposed that a quantum analog of liquid crystal states - electronic smectic and nematic 
phases - may play an important role in the physics of strongly correlated systems [76] , such as high temperature 
superconductors. The smectic and nematic phases are characterized by the broken unidirectional translational and 
rotational symmetry and broken rotational symmetry, respectively. In cuprates, the electronic liquid crystal (ELC) 
phases may have different origins. Starting from a hole doped AF Mott insulator, electronic smectic or nematic 
phase could correspond to static or fluctuating stripes [13,14]. They can also correspond to static or fluctuating 
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Figure 11. Cartoon picture suggesting the indication of a spin and charge stripe structure and an orbital-like magnetic order. The 
stripe structure is characterized by doped holes located on two-leg ladders, separating AF domains. The charge periodicity is ~ 6a 
and the magnetic periodicity is ~ 12a. 

spin spiral phases [77]. Starting from the metalhc side of cuprates phase diagram, an electronic nematic phase can 
take the form of a spontaneous distortion of the Fermi surface [78], breaking the a-b symmetry. 

Direct or indirect evidence of electronic smectic and nematic phases have been reported in superconducting 
cuprates, via different experimental techniques [76]. For instance, in strongly underdoped YBC06.45 [19], quasi- 
ID incommensurate (IC) spin fluctuations at Q/c=(0.5=b5, 0.5) spontaneously develop below ^ 150 K and become 
static at very low temperature. The observation can be generalized for hole doping below 9% (corresponding to 
an oxygen content < 0.5 in YBCOe+x) [20]. In strongly underdoped YBCO, both neutron scattering [19,20] and 
resistivity [79] measurements indicate that the a-b symmetry is broken. These observations strongly support the 
existence of ELC phases in strongly underdoped YBCO, but there is not yet a consensus concerning the electronic 
instability at the origin of the phenomenon. In twinned LSCO samples, IC spin fluctuations at Qic = (0.5di^, 0.5) = 
(0.5, 0.5 ± 5) are also observed and are usually interpreted in terms of fluctuations stripes oriented along Cu-0 
bonds (either along a or b due to the existence of twin domains). The nematic ELC phase is therefore associated 
with the appearance of these fluctuating stripes, breaking the a-b symmetry of the system [13,18]. Obviously, 
these ELC phases and the novel Q=0 magnetic order correspond to different electronic instabilities. This leads a 
simple question: what happens when these phases met in the phase diagram of cuprates ? 

To address this issue, it is first interesting to consider the case of LSCO where we have reported a short range 
Q=0 AFO[36]. Interestingly, low energy IC spin fluctuations display a feedback to the short range ordering. 
The incommensurability parameter 6 is renormalized at Tmag- The original picture of stripes physics implies 
that the doped holes self-organise into lines of charged stripes which create antiphase domain walls between 
antiferromagnetic domains where Cu^^ spins fluctuate and freeze at low enough temperature. A sketch of such 
charge and spin organisation is given in Fig. 11 for a hole doping of 0.085 ^ 1/12 corresponding to our LSCO 
sample [36]. We here consider the case of bond-centered stripes where a doped hole is located on two-leg ladders. 
The orbital-like order could be confined within the two-leg ladders, yielding typical correlation lengths of ^ 2a 
along the (100) direction and ^ 3a along the (110) direction, as observed experimentally (see Figs. 9. A and B). 
The orbital-like magnetic order should not develop a phase coherent from one ladder to another, otherwise the 
magnetic intensity should appear at incommensurate wave vectors away from H=l, instead of the broad peak 
observed around Q=(1,0,L). This example in LSCO underlines the competition between ELC and orbital-like 
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Q=0 instabilities. 

In YBCO, neutron scattering measurements provides evidence of ELC phases, but only in the strongly under- 
doped regime (p< 0.09). The Q=0 AFO has been studied from YBCOe.s (p=0.09) up to (Y,Ca)BC07 (p=0.2). 
In this hole doping range, neutron scattering studies do not indicate that the ELC phases dominate the low 
energy properties of the system. In YBCOe.e (p=0.12) for instance, the spin excitation spectrum exhibits a spin 
pseudogap in the normal state and low energy spin fluctuations are hardly measurable. Above the spin pseudogap, 
IC spin fluctuations still display an a-b anisotropy in the energy range between ~20 meV and ~38 meV [80]. To 
date, one still ignore at which temperature this anisotropy appears. In the hole doping range where the Q=0 
AFO has been identified, an a-b anisotropy has been recently detected in the Nernst effect [81]. In principle, this 
supports theoretical approaches [13,14] in which the C4 symmetry is spontaneously broken in the pseudogap state, 
owing to strong electronic correlations yielding the nematic ELC state. This anisotropy develops when entering 
the pseudogap state at a temperature corresponding exactly to Tmag of the Q=0 order in a large doping range. 
Clearly, the anisotropy starts at a temperature where there is no indication of a similar anisotropy neither in the 
spin excitation spectrum nor in electrical resistivity. In addition, the existence if an a-b anisotropy in the Nernst 
effect has not been studied in strongly underdoped YBCO, where the ELC phases take place. Actually, it exists 
a completely different interpretation of this anisotropic Nernst effect [82]. In this interpretation, the anisotropic 
Nernst effect would appear as a consequence of the time-reversal symmetry breaking below T* in agreement with 
polarized neutron data [33,34] and Kerr effect [46]. 

In underdoped YBCO for p>0.09, quantum oscillations data indicate the existence of small electron Fermi 
pockets [22,23,24] at very large magnetic field (H>40 Tesla) when superconductivity is supposed to be destroyed. 
These data are often interpreted in terms stripe- like spin density wave [25]. Recent transport [83] and quantum 
oscillations [84] data further indicate that the existence of these electron Fermi pockets is quite sharply peaked 
around the hole doping of 1/8, where stripes are expected to be pinned down by the lattice. One can be tempted 
to establish a connection between the ELC phase reported in strongly underdoped YBCO (for p<0.09) and the 
small Fermi pockets reported at larger doping (p> 0.09) and very high magnetic field (H>40 T). This interplay is 
further supported by the enhancement of the smectic ELC phase in YBCO6.45 when applying a modest magnetic 
field of 14 Tesla [21]. Strictly speaking, however, there is no compelling evidence of well developed ELC phases in 
the hole doping range where the quantum oscillations are reported, questioning this straightforward scenario. 

At variance with LSCO, in YBCO there is no clear indication that the Q=0 AFO and ELC phases could be 
found in the same portions of the phase diagram versus hole doping and magnetic field. In order to get better a 
understanding of the interplay between these phases, the search for the Q=0 magnetic order should be extended to 
the strongly underdoped regime. The search for such an order at high magnetic field is for the moment technically 
impossible with polarized neutron diffraction technique. However, instead of applying a uniform perturbation to 
the system, one should consider the effect of local perturbations induced by the substitution of non magnetic Zn 
or magnetic Ni impurities. If there were a competition between the Q=0 magnetic order and the ELC phase, their 
balance could be modified around local defects. 



5. Conclusion 

The elucidation of the pseudogap phenomenon of the cuprates has been a major challenge in condensed matter 
physics for the past two decades. Following initial indications of broken time-reversal symmetry in photoemission 
experiments [29], recent polarized neutron diffraction works demonstrated the universal existence of an unusual 
magnetic order below T* in 3 different cuprates families [33,34,35,36] whatever the number of Cu02 plane per unit 
cell. The observed magnetic order can be interpreted as a Q=0 antiferromagnetic order preserving the translational 
symmetry of the lattice. The observed symmetry is furthermore consistent with a particular type of order involving 
circulating orbital currents (CC — On phase) [10,11], and with the notion that the phase diagram is controlled by 
a quantum critical point [66]. The ordered temperature indeed extrapolates to zero at the doping level of p ~ 0.19 
corresponding to the end point if the pseudogap state. To date, it is the first direct evidence of a hidden order 
parameter associated with the pseudogap phase of high-Tc cuprates. It has the profound implication that the 
pseudogap transition T* is a genuine phase transition although no sharp thermodynamic anomalies are observed 
[45]. Uniform susceptibility [48] and magneto-optic Kerr effect [46] measurements confirm this picture. Recent 
measurements in La2-xSrxCu04 reveal that a similar ordering occurs in that system for x=0.085 [36]. However, 
the observed order is at short range and bidimensional. That difference might be related to the competition of 
the Q=0 magnetic order with the stronger tendency towards an electronic liquid crystal state in that system. 
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